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Abstract 

It was observed by Tod [2] and later by Dunajski and Tod [I] that the Boyer-Finley 
(BF) and the dispersionless Kadomtsev-Petviashvili (dKP) equations possess solutions whose 
level surfaces are central quadrics in the space of independent variables (the so-called central 
quadric ansatz). It was demonstrated that generic solutions of this type are described by 
Painleve equations Pjjj and P/j, respectively. The aim of our paper is threefold: 

- Based on the method of hydrodynamic reductions, we classify integrable models possessing 
the central quadric ansatz. This leads to the five canonical forms (including BF and dKP). 

- Applying the central quadric ansatz to each of the five canonical forms, we obtain all 
Painleve equations Pj — Pyx , with Pyi corresponding to the generic case of our classification. 

- We argue that solutions coming from the central quadric ansatz constitute a subclass of 
two-phase solutions provided by the method of hydrodynamic reductions. 

MSC: 34M55, 35C05, 35L10, 35Q51, 35Q75, 37K10. 

Keywords: multidimensional integrability, hydrodynamic reductions, central quadric ansatz, 
Painleve equations. 
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1 Introduction 



This paper grew from an attempt to understand the construction of [144 1] reducing the Boyer- 
Finley (BF) and the dispersionless Kadomtsev-Petviashvili (dKP) equations to Painleve trans- 
cendents via the so-called 'central quadric ansatz'. This procedure applies to PDEs of the form 

(a(u)) xx + (b(u)) yy + (c(u)) tt + 2(p(u)) xy + 2(q(u)) xt + 2{r{u)) yt = 0, (1) 

and consists of seeking solutions u{x, y, t) in implicit form, 

(x,y,t)M{u){x,y,t) T = 1, (2) 

where M(u) is a 3 x 3 symmetric matrix of u. The level surfaces of such solutions, u = const, are 
central quadrics in the space of independent variables x, y, t. Remarkably, for both BF and dKP 
equations considered in [Hid], differential equations for the matrix M{u) reduce to Painleve 
transcendents. Although the occurrence of Painleve transcendents via similarity reductions 
of soliton equations is a well-known fact, their appearance in the context of multidimensional 
dispersionless integrable PDEs is, to the best of our knowledge, an entirely new phenomenon. 
In this paper we address the following questions: 

• Classify integrable PDEs of the form ([1]). Our approach to this problem is based on the 
method of hydrodynamic reductions which, as demonstrated in [1], provides an efficient 
criterion for the classification of multidimensional dispersionless integrable equations. In 
particular, both BF and dKP equations are known to possess an infinity of hydrodynamic 
reductions. The classification is to be performed modulo (complex) linear changes of the 
independent variables x,y,t, as well as transformations u — > f(u), which constitute the 
equivalence group of our problem. 

• Describe Painleve reductions of integrable PDEs resulting from the above classification by 
applying the central quadric ansatz. 

• Relate the central quadric ansatz to the method of hydrodynamic reductions. 
To formulate our first result we introduce the symmetric matrix 

(a' p' q' 

p' b' r' 

q' r' c' 

where prime denotes differentiation by u. 

Theorem 1 A PDE of the form {IJJ is integrable by the method of hydrodynamic reductions if 
and only if the matrix V{u) satisfies the constraint 

V" = (lndetV)V + kV, (3) 

for some scalar function k. Modulo equivalence transformations, this leads to the five canonical 
forms of nonlinear integrable models possessing the central quadric ansatz: 

UXX + Uyy ~ g{u)yy ~ ^ (u) U = 0, 
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U XX + Uyy + (e U ) tt = 0, 

(e M - u) xx + 2u xy + (e M )« = 0, 

(^ 2 )xy + %j/ + 2Ua;t = 0, 

/iere g(u) = ln(l — e u ). Examples 2 and 4 are the familiar BF and dKP equations. 

We point out that the constraint ©, which implies V" G span{V, V'}, means that the 
'curve' V(u) lies in a two-dimensional linear subspace of the space of 3 x 3 symmetric matrices. 
The classification of normal forms of such linear subspaces leads to the five canonical forms of 
Theorem 1. Equations listed in Theorem 1 are not new: in different guises, they have appeared 
in the classification of multidimensional integrable systems (see Sect. 2 for the proof of Theorem 
1 and further discussion). 

In Sect. 3 we apply the central quadric ansatz ([2]) to the five canonical forms of Theorem 
1. Differentiating ([2]) implicitly and substituting into ([T]) one obtains equations for the matrix 
M(u) which can be most naturally written in terms of the inverse matrix N = — M _1 : 

gN' = V where g 2 det A £ const, (4) 

see |144 [T]. For matrices V associated with the five canonical forms of Theorem 1, generic 
solutions of equations @) reduce to the following Painleve transcendents: 

Case 1: General case of Pyi', 

Case 2: Special case of Py reducible to Pm; 

Case 3: General case of Py; 

Case 4-' General case of Pjj with a reduction to Pi; 
Case 5: General case of Pjy. 

It is quite remarkable that the whole variety of Painleve equations is contained in the pair 
of matrix equations , @) : setting £ = 1 and eliminating the scalar g we can rewrite them in 
the equivalent form, 

N' = VVdet N, (V/ det V)' = kV. 

In Sect. 4 we show that, at least for the dKP equation, solutions coming from the central 
quadric ansatz constitute a particular subclass of two-phase solutions provided by the method 
of hydrodynamic reductions. 



2 Classification of integrable PDEs of the form (CD): proof of 
Theorem 1 

Following the method of hydrodynamic reductions [H H], we seek multi-phase solutions of 
equation ([1]) in the form 

u = u(R 1 ,...,R N ), (5) 
where the phases R l (x, y, t) are required to satisfy a pair of compatible systems of hydrodynamic 

type, 

Ri = X i (R)Rl Ri = ^(R)Ri, (6) 
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R = (R l , . . . , R N ). We say that equations ([U]) provide an iV-component hydro dynamic reduction 
of the equation ([1]). Note that the number of phases, N, is allowed to be arbitrary. The 
compatibility of equations ([6]) requires the following conditions [15], 



A* _ /4 



(7) 



i j, here A*- = d^X 1 , ju*- = d R jfi l . It was observed in [3] that the requirement of the existence 
of iV-phase solutions parametrized by N arbitrary functions of one variable imposes strong 
constraints on the original system ([1]), and provides an efficient approach to the classification of 
multidimensional dispersionless integrable systems. 

In the present context the integrability conditions can be derived as follows. Upon substitu- 
tion of © and ([6]) into equation ([I]) we obtain that the characteristic speeds A* and fji 1 obey the 
dispersion relation, 

A' 4 = a' + b'fi i2 + c'X l2 + 2 (p'fj + q'X' 1 + r'AV ) = 0. (8) 

Differentiating this relation with respect to B? , i ^ j, and taking into account the equations © 
and (JSJ), we obtain 

V = lW-*)§j^ 4 = ltf-rt£j"i> (9) 
where we introduced the notation 

A* = a" + 6'V 2 + c"X i2 + 2 (pV + ^''A* + r"AV*) , 

A ij = a' + feVV J ' + c'X l X j + p'(Li l + Li j ) + q'(X i + X j ) + r'(AV + A J V)- 

Moreover, we obtain the following expressions for the second order derivatives of u, 

A* 

where 



UiUj, (10) 



A^' = a" + &W + ^'AU 3 ' + p'V + + g"(A* + A J ) + r"(AV + *V 



Equations (J9j) and (|10p constitute the so-called generalised Gibbons-Tsarev system governing 
multi-phase solutions (hydrodynamic reductions) of equation ([T|). The conditions of compat- 
ibility of Eqs ©-CEO]), (Xj)k = (X\)j, (fj,j)k = (fJ.\)j and (uij) k = (uik)j, lead to polynomial 
expressions in the characteristic speeds A\A J , A fc and /x 4 ,^ , fi k which are required to vanish 
modulo the dispersion relation ([8|). This leads to the integrability conditions ([3]). 

To solve Eq. ([3]) let us view V(u) as a curve in the space of 3 x 3 symmetric matrices. Since 
the acceleration vector V" is a linear combination of the position vector V and the velocity 
vector V', this curve must belong to a fixed two-dimensional linear subspace (the case when this 
subspace is one-dimensional corresponds to linear equations: in this case the central quadric 
ansatz was analysed in [2]). The classification of two-dimensional linear subspaces of the space 
of 3 x 3 symmetric matrices leads to the five canonical forms which are presented below in the 
format span{A, B} for some particular choices of symmetric matrices A and B. Recall that pairs 
of symmetric matrices, A and B, are classified by Jordan normal forms of the operator BA^ 1 , 
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further subcases correspond to coincidences among eigenvalues of this operator. Choosing a 
suitable basis A, B in the subspace span{A, B} leads to the five canonical forms mentioned above. 
In each of these cases one can set V = A + f(u)B (to be precise, V = g(u)A + f(u)B, however, 
g(u) can be set equal to 1 via an equivalence transformation u — > (f{u)). The substitution into 
the integrability condition ([3]) leads to simple ODEs for f(u) which can be readily solved. 

Case 1. 

A 

The substitution of V = A + f(u)B into the integrability condition ([3]) leads to k = and 
/" = f' 2 (^j + JTT^ ' so that without any loss of generality one can take f(u) = e u /(l — e u ). 
This leads to the first equation of Theorem 1. Setting u = v yy + vu we obtain the equation 
e v X x+v yy _|_ e v t t+vyy _ ^ QJ . ^ e q U i va lently, e Vxx + e Vtt = e~ Vyv , which first appeared in [5] in 
the classification of integrable PDEs of the form F(v xx ,v yy ,vu) = 0. It can be viewed as an 
integrable generalisation of the BF equation. 

Case 2. 
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The substitution of V = A + f(u)B into the integrability condition ([3]) leads to k = and 
/" = /' 2 //i so that without any loss of generality one can take f(u) = e u . This leads to the BF 
equation. 

Case 3. 

A 

The substitution of V = A + f(u)B into the integrability condition ([3]) leads to k = and 
/" = Z' 2 / (/ + !)> so that without any loss of generality one can take f(u) = e u — 1. This leads to 
the third equation of Theorem 1 which is yet another integrable deformation of the BF equation: 
its equivalent forms, u xy = (e u )u — c(e u ) xx and u xy = (e u )u + (e u )t y , have appeared in [HE]. Both 
can be reduced to case 3 via appropriate (complex) linear transformations of the independent 
variables x, y, t. An alternative first order form of this equation is provided by the system 
v t = Vy v +Z V » Wx = VV v+w w hich appeared in [B] in the classification of integrable Hamiltonian 
systems of hydrodynamic type in 2+1 dimensions. Setting v = s xy , w = st y we obtain the second 
order PDE s xy + st y = e Sxt . For u = s x t this gives the equation u xy + ut y = {e u ) x t which is also 
equivalent to case 3. 

Case 4. 

A 

The substitution of V = A + f{u)B into the integrability condition ([3]) leads to k = and 
/" = 0, so that without any loss of generality one can take f(u) = u. This leads to the equation 
(uu x ) x + 2u xy + uu = 0, which is equivalent to the dKP equation. 
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Case 5. 
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The substitution of V = A + f(u)B into the integrability condition © leads to = and 
/" = 0, so that without any loss of generality one can take f(u) = u. This leads to the last 
equation in Theorem 1. Setting u = v xy we obtain the equation v xy + v yy + 2v x t = 0. In this 
form, it appeared in [13J as the simplest case in the classification of integrable hydrodynamic 
chains satisfying the so-called 'Egorov' property. 
This finishes the proof of Theorem 1. 



3 Reduction to Painleve equations 

In this section we investigate systems ©, 

gN' = V, g 2 det N = f = const, 

for all canonical forms of Theorem 1, and reduce them to Painleve transcendents. For the BF 
and dKP equations this was done in |14|H]: we include these results for the sake of completeness. 
We find it more convenient to work with a- forms of Painleve equations as listed in [11] , Appendix 
C. 

Case 1: u xx + u yy — ln(l — e u ) yy — ln(l — e u )u = 0. We have 
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so that the first equation (jj]) implies that a, j3, 7 are constants, while F' = |, G' = H' = 

j^y g . Thus, G' = j^F', H' = j^F>. Setting s = ^ one obtains G s = sF s , H s = 
(s — l)F s . Parametrizing these relations in the form F = o~ s , G = so~ s — o~,H = sa s — a + \i — a s , 
and substituting into the second equation 

FGN - a 2 H - (3 2 G - ^ 2 F + 2a/3 7 = £F' 2 , 

(note that F' = s(s — 1)F S ), we obtain an ODE for o~(s), 

£s 2 (s - l) 2 a 2 ss = a s (sa s - a) (sa s -a + fi-a s )- ^ 2 a s - j3 2 (sa s -a) -a 2 (sa s - a + /i - o s ) + 2a/3j. 

This ODE belongs to the class 

s 2 (s - l) 2 a 2 ss = aa s (sa s - a) 2 + (sa s - a)P 2 (a s ) + Pi(a s ), 

where Pi and Pi are polynomials in a s of the order two and one, respectively, such that the 
leading term of P2 is —aa 2 . Modulo transformations a — > a\o + a 2 this form is equivalent to 
the general case of Pvi, 

a s s 2 (s - l) 2 a 2 s + 2[a s (sa s - a) - a 2 - v\v 2 v z v$ = (a s + vf)(a s + i/f ) (<r s + vf){(T s + v\). 

see m. 
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Case 2: u xx + u yy + (e u )tt = 0. We have 
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so that the first equation Q implies that a, /3,7 are constants, while F' = l/g, G' = 1/g, H' = 
e u /g. Thus, G' = F', F' = e u F'. Setting e u = s, parametrizing these relations in the form 
F = a s , G = a s + p, H = sa s — a, and substituting into the second equation @ we obtain an 
ODE for a(s), 

£sV 2 s = a s (a s + p)(sa s - a) - j 2 a s - a 2 {so s - a) - f3 2 (a s + pi) + 2a(3j. 
This ODE belongs to the class 

s2<j2 ss = ( S(7 s ~ a)P 2 (a s ) + Fi(ct s ), 

where P 2 and P\ are polynomials in <r s of the order two and one, respectively. This is the special 
case of Py which is reducible to Pni, see Case 3 below. 

Case 3: (e u - u) xx + 2u xy + (e u ) u = 0. We have 



/ e n - 1 1 \ 
V = 1 0,1= 
\ e" / 
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so that the first equation @ implies that a, /3, 7 are constants, while F' = (e w — l)/<?, G' = 
1/c/, JET' = e u /g. Thus, F' = e u G', F' = (e u - 1)G'. Setting e u = s, parametrizing these 
relations in the form G = <r s , F = sa s — a, F = sa s — a + p, — a s , and substituting into the 
second equation Q we obtain an ODE for <r(s), 

£sV 2 s = /3(scr s — a + pt — a s )(sa s - a) - j 2 (sa s - a + p, — a s ) — a 2 (sa s - a) + 2aja s - a 2 /3. 

This ODE belongs to the class 

s 2 a 2 s = a(sa s - a) 2 + (sa s - a)P 2 (a s ) + Pi(a s ), (11) 

where F2 and P% are polynomials in a s of the order two and one, respectively. In the case 
transformations a — > a\o + a 2 s + 03, s — > b\s bring any generic equation of the form (|11|) to the 
canonical Py form, 

s 2 cj 2 s = [u- sa s + 2a 2 + (yx + v 2 + i>z)<Js\ 2 ~ 4cr s (i/i + <r s ){v 2 + a s )(u 3 + a s ). 
The case a = leads to the special case of Py which is reducible to Pni, 

s 2 a 2 ss = (sa s - a)P 2 (a s ) + Pi{a s ). 



Case 4: u x t — (uu x ) x — u yy = 0. We have 

/ -u 1/2 \ (Fa G \ 

F= -1 , iV = a F /3 , 

V 1/2 J \ G /3 7 / 
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so that the first equation © implies that a,/?, 7 are constants, while F' = —u/g, G = 
l/(2g), H' = —1/g. Thus, G' = —H'/2, F' = uH' . Parametrizing these relations in the 
form H = a', H = \i — a' /2, F = ua' — a, and substituting into the second equation @ we 
obtain an ODE for a(u), 

£( a ") 2 = ja'(ua' -a)- /3 2 (ua' -a)- a'(u - a'/2) 2 + 2aj3(u - a'/2) - a 2 7 . 

This ODE belongs to the class 

(a") 2 = aa'(ua' - a) + b(ua' - a) + P 3 (<A (12) 

where P3 is a third order polynomial in a' . In the case one can eliminate b by a transfor- 

mation a — > a + su. Further transformations a — > a±a + a>2, u — > b\u + 62 bring any equation 
of the form (|12p to the canonical Pjj form, 

(a") 2 + A(a'f + 2a' {ua' - a) + 8 = 0. 

Similarly, for a = 0, 67^ 0, transformations a — > a\a + su + a^, u — > &1U + &2 bring any equation 
of the form (|12p to the canonical Pj form, 

(a") 2 + 4(a'f + 2(ua' - a) = 0. 
Case 5: (u 2 ) xy + u yy + 2u x t = 0. We have 
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so that the first equation (J3|) implies that a,/3, 7 are constants, while G' = 1/g, H' = 1/g, F' = 
u/g. Thus, H' = G' , F' = uG' . Parametrizing these relations in the form G = a' , H = 
a' + /j,, F = ua' — a, and substituting into the second equation @ we obtain an ODE for a(u), 

£(a") 2 = - 7 (W - a) 2 + 2j3a' (ua' -a)- a' 2 (a' + n) + aj(a' + fi) - a/3 2 . 

This ODE belongs to the class 

(a") 2 = a{ua' - a) 2 + ba '(ua' - a) + P 3 (a'), (13) 

where P3 is a third order polynomial in a'. In the case a^0 one can eliminate b by a translation 
of u. Further transformations a — > a\a + a2U, u — > b\u bring any equation of the form ()13|) to 
the canonical Pjv form, 

(a") 2 = A(ua' - a) 2 - 4a' (a' + ^)(a' + u 2 ), 

Similarly, for a = 0, b 7^ 0, transformations a — > a\a + 03, u — > b\u + b% bring any equation of 
the form (|13p to the canonical Pu form, 

(a") 2 + 4(a'f + 2a' (ua' -a) +8 = 0. 
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4 Central quadric ansatz versus hydrodynamic reductions 

The main observation of this section is that solutions coming from the central quadric ansatz 
([2]) can also be obtained as a particular case of two-phase solutions, 



u 



u(R l ,R 2 ), 



where the phases R 1 (x, y, t), R 2 (x, y, t) satisfy a pair of two-component systems of hydrodynamic 
type, 

Rl = X\R)R l x , Ri = fj,\R)Ri, (14) 

i = 1,2. The general solution of equations (fT4"j) is given by the so-called generalised hodograph 
formula [15], 

x + \ 1 (R 1 ,R 2 )t + u 1 (R 1 ,R 2 )y = v 1 (R 1 ,R 2 ), 

x + X 2 (R\ R 2 )t + u 2 (R\R 2 )y = v 2 (R\R 2 ), [ j 

where v 1 are the characteristic speeds of yet another flow, 

K = u\R)Rl (16) 

which commutes with (fTl"]) . Equations (fT5"j) define R 1 (x,y,t) and R 2 (x,y,t) implicitly. These 
relations can also be viewed as the equations of a line congruence (two-parameter family of 
lines) in the space of independent variables x, y, t, parametrised by R , R 2 . Asw = u{R 1 , R 2 ) is 
constant along the lines of this congruence, the level surfaces u = const will automatically be 
ruled. Since a quadric carries two (complex) one-parameter families of lines, solutions constant 
on quadrics can therefore be viewed as two-phase solutions in 'two different ways': they are 
constant along two distinct congruences of lines. 

Taking dKP equation as an example, we will now provide explicit derivation of the central 
quadric ansatz from the general two-phase solutions. Let us rewrite dKP as a two-component 
first order system, 

U t = UU X + Wy, Uy = W X . 

It will be more convenient for our purposes to work with variables u, w rather than R 1 , R 2 . 

Proposition 1. Two-phase solutions of the dKP equation, u(x,y,t) and w(x,y,t), are given by 
the implicit formulae 

x + (z u + u)t = m u , y + z w t = m w , (17) 
where the functions z(u, w) and m(u, w) satisfy the PDEs 

Zuu z w z uw z u z ww + 1 — 0, (18) 



Proof: 

The analogues of equations (|14p read 

u t = uu x + z x , w t = q x and u y = w x , w y = z x , (19) 

where the commutativity conditions, ut y = u y t and wt y = w y t, imply that the functions q(u,w) 
and z(u,w) satisfy the relations q u = z u z w , q w = z^ + z u + u. Notice that the compatibility 
condition, q uw = q wu , implies the first relation (fT8|) . This equation for z{u,w) first appeared in 
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[7j in the context of hydrodynamic reductions of Benney's moment equations. It is known as 
the Monge- Ampere form of the Gibbons-Tsarev system. The analogue of the third commuting 
flow (USD is 

u T = m x , w T = n x , (20) 

where the commutativity conditions imply that the functions m(u, w) and n(u, w) satisfy the 
relations n u = z u m w , n w = z w m w + m u . Their compatibility condition, n uw = n wu , implies the 
second relation (|18p . Rewriting equations (|19|) and (|20p in matrix form we obtain 

u \ ( z u + u z w \ ( u \ ( u \ ( 1 \ f u 



w J t \ z u z w z 2 w + z u + u J \ w J x \w J \ .:„ J V 

and 

u \ = ( m u m w \ ( u 

w ) r \ z u m w z w m w + m u J \ w 

respectively. In matrix notation, the generalised hodograph formula (|15p is 

l0\ + fz u + u z w \ t+ f° 1 \ = ( m u m w 

1/ y ^u^w %w ~i~ ^ / V %u Z w J y Z u ni w Z w 71l w -\- 71 1 u 



One can readily see that only two of these relations are independent, and they are equivalent to 

Remark. Differentiating ()17p implicitly by x,y,t we obtain the relations 

u t = (z u + u)u x + z w u y , w t = (z u + u)w x + z w w y , (21) 

which manifest the invariance of the corresponding solutions under the conditional symmetry 

X = (z u + u)d x + z w d y - d t . (22) 

The presence of such symmetry is characteristic of two-phase solutions for equations in 2 + 1 
dimensions [9], [10] ■ Differentiating (|2ip with respect to x,y,t and eliminating the derivatives of 
z we obtain second order differential constraints for u and w characterising two-phase solutions, 

&xy u tt + $ x t u yy + 8 yt u xx - 2 (5 xt 5 X yU t y + 5 tx 5tyU xy + 

) = o> 

SlyWU + Sl t Wyy + XX ~ 2 {5 xt 8 X yW t y + S tx 5 t yW X y + 5y t 5y X W tx ) = 0, 

where the functions 5 y t,8t x ,5 xy are 2x2 minors of the Jacobi matrix J = d(u,w)/d(x,y,t), 
i.e. 5 y t = UyWt — utWy, etc. Introducing the vector field X = (5 y t,5t x ,5 xy ) one can rewrite 
these constraints in the equivalent form XUX T = 0, XWX T = where U, W are the Hessian 
matrices of u and w. The constraints (|23|) completely characterize two-phase solutions of the 
dKP equation, providing an easy criterion to verify if a given solution can be constructed using 
this method. We emphasize that differential constraints (|23p are universal: they govern two- 
phase solutions of any two-component quasilinear system. Indeed, the general solution of (f23|) 
is given by implicit relations 

x + p(u,w)t = q(u,w), y + r(u,w)t = s(u,w), (24) 
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where p, q, r, s are arbitrary functions. In other words, functions u, w solve (|23|) if and only if they 
are constant along a two-parameter family of lines. This follows from the following equivalent 
geometric representation of the constraints ([23]): 

L x u = L x w = 0, XUX T = XWX T = 0. 

The first two conditions mean that the functions u and w are constant along integral trajectories 
of the vector field X. This, in particular, implies that X = (8yt,St x ,S xy ). The last two conditions, 
which are equivalent to (j23|) . mean that integral trajectories of X are asymptotic curves on the 
level surfaces u = const and w = const. It remains to point out that if two surfaces intersect 
transversally along a curve which is an asymptotic curve on both of them, this curve must be a 
straight line (indeed, its osculating plane is tangential to both surfaces, and hence degenerates 
into a line). Adding the relations (|24|) to the dKP equation one obtains constraints for p,q,r,s 
resulting in (fT7|) . 

Following [T2] we will demonstrate that the system (JTHj) possesses solutions expressible in 
terms of the Painleve equation Pjj. 

Proposition 2. The system ilty) possesses solutions of the form 

z(u, w) = si(u) + S2{u)e w + s 3 {u)e~ w ', m(u, w) = S2(u)e w — s 3 (u)e~ w , 

where the functions si(u), 82(11), s 3 (u) satisfy the ODEs 

s[ = 2s 2 s 3 - u, s 2 ' = (2s 2 s 3 - u)s2, s'g = (2s 2 s 3 - u)s 3 . (25) 

This expression for z was obtained in [12J based on the method of linear differential constraints. 
Substituting the above z, m into the equations (|19j) and eliminating w, one obtains an implicit 
relation of the form ((2J) for u, 

{s 2 s 3 )' 2 {t 2 - 1) - s 2 s 3 (x + 2s 2 s 3 tf + (s 3 s' 2 - s 2 s' 3 ){x + 2s 2 s 3 t)y + s' 2 s' 3 y 2 = 0, (26) 

which shows that u is indeed constant on central quadrics. It remains to demonstrate that the 
system ([25]) is equivalent to Pjj. Indeed, equations (j25|) possess a conservation law, s 3 s' 2 — S2s' 3 = 
c = const, which can be resolved in the form s' 2 = ps2 + c/s 3 , s' 3 = ps 3 . This, in particular, 
implies the relation (s2S 3 )' = 2ps2S 3 + c. Differentiating s 3 = ps 3 and using the last equation 
(|25|) we obtain 2s2S 3 = p'+p 2 +u. Substituting this expression for S2S 3 into the previous relation 
we obtain the Pjj equation, 

p" = 2p 3 + 2pu + 2c - 1. 
Ultimately, all coefficients of the relation (|26p are expressed in terms of p, 

S2S 3 = ~{p' + p 2 + u), s 3 s' 2 - s 2 s 3 = c, s' 2 s 3 = p 2 s 2 s 3 + cp= ^p 2 (p +p 2 + u) + cp. 
This gives a solution equivalent to the one constructed in pQ. 
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